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Abstract. A Brownian information machine extracts work from a heat bath through 

a feedback process that exploits the information acquired in a measurement. For the 
O ' paradigmatic case of a particle trapped in a harmonic potential, we determine how 

r^ , power and efficiency for two variants of such a machine operating cyclically depend 

on the cycle time and the precision of the positional measurements. Controlling only 
2 ' the center of the trap leads to a machine that has zero efficiency at maximum power 

^ I whereas additional optimal control of the stiffness of the trap leads to an efficiency 

bounded between 1/2. which holds for maximum power, and 1 reached even for finite 

cycle time in the limit of perfect measurements. 
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In Kelvin's formulation, the second law of thermodynamics states that no work 
can be extracted from a thermally equilibrated system through a cyclic process that 
leaves no trace elsewhere. If, however, more detailed information of the system becomes 
available through a measurement, then one can indeed extract work as illustrated a long 
time ago with the gedankenexperiments of Maxwell's demon and Szilard's engine [1]. 
More recently, by combining concepts from stochastic thermodynamics with those from 
information theory, a quantitative framework has emerged leading to bounds refining 
the second law to such feedback driven processes [SKSHUlSlEllTllHlEl [in]. Specialized 
to one cyclic process starting in equilibrium, the bound 

W <X (1) 

connects the mean extractable work W to the mean information I (defined more 
precisely below) acquired through the measurement. Brownian particles in time- 
dependent potentials provide a paradigm for such systems both in recent experiments 
[m [12] as in several theoretical case studies [131 El CHj- The latter works have 
demonstrated that saturating the bound in ([T|) typically requires both an infinite cycle 
time and a sufficient number of control parameters in the potential. 

The purpose of the present paper is to study these processes from a perspective that 
focuses on the performance of such Brownian information machines in a steady state 
where measurements and subsequent optimal driving based on these are repeated with a 
finite cycle time t. On average per cycle, by exploiting the information X,,, the machine 
extracts the work VT* thus delivering a power P = W.^/t. The extant generalization of 
the bound ([T|) to such a cyclic operation [8] then motivates to define efficiency as 

7] = W/I (2) 

following in spirit an earlier approach [1]. Apart from maximum efficiency, it is 
particularly interesting to determine efficiency at maximum power. The later concept 
has been studied extensively for non-feedback driven heat engines operating between 
two heat baths, see, e.g., [Ml [T71 [181 [19] and references therein, and, more recently, also 
for autonomous isothermal machines |20j . 

The solution to this problem of performance at a finite cycle time cannot trivially 
be inferred from available results [H] on the maximal extractable work following one 
measurement in finite time since at the beginning of the second (and any further) cycle 
the system will typically not have reached thermal equilibrium again. In fact, the initial 
state of the z— th-cycle will depend on the result of all previous measurements which 
makes the present problem non-trivial. 

Our system consists of an overdamped Brownian particle in a harmonic potential 

V{x,T) = k{T)[x-\{T)]y2 (3) 

with external time- dependent control of the center, A(r), and stiffness, A;(r), of the trap 
[T^ . Throughout the paper, we use dimensionless variables. The harmonic potential 
has the advantage that a Gaussian distribution 

PW=^.-(M^)-p^exp(-(H^) (4) 
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remains Gaussian both under the stochastic dynamics in the potential and under 
positional measurements with an error ± y^- The dynamics of the mean fe(r) and 
variance ^^(t) of x follows from the corresponding Fokker-Planck equation as [11] 

h{r) = k{r)[\{r)-h{r)] (5) 

and 

y{r)=y{r)[lly\r)-k{r)] (6) 

where we denote time-derivatives with a dot throughout. 

We now implement a cyclic feedback scheme based on measurements of the position 
repeated periodically in intervals of lengths t. At the beginning of the i— th cycle, we 
measure the position Xj with a precision ± y^ leading to the distribution 

p{X,)=UxM.{yif + yl:) (7) 

for the measured value if the distribution prior to the measurement is characterized by 

p-(x)=A4(6r, (%)')• (8) 

After the measurement the distribution for x follows from Bayes' theorem as 

p{x\x,)=UMAytf). (9) 



with 



and 



ht = ''^^^\+^y- (10) 

ivi ) + y?, 



^m 



(y+f = ^y^] y'r- . (11) 

{Vi ) + yli 

Based on this measurement, we maximize the extracted work by optimally adjusting 
the control parameters. Quite generally, given an initial state 6(0) = b^,y'^{0) = {yf) 
and a time-dependent hij) and ?/(r), the extracted work after a time t becomes [H] 

with 

- Wr = [h\t) - {btf]/2 + fdr h\r) (13) 



^0 

and 

- 1^7* = [y\t) - {ytf]/2 - Hy{t)/yt] + f dr fir). (14) 

Here, we have required that the trap is centered at A = with stiffness k = 1 aX 
beginning and end of the cycle allowing for jumps of these two control parameters. 
Depending on the amount of control available, two cases must be distinguished. 

If the stiffness is fixed, k^r) = 1, only the center of the trap A(r) is controllable. 
Using (Q in the integral (TT4l) shows that in this case W°^^ = 0. The fe-dependent term is 
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maximized by a linear function 6(r) = ^^^[1 — r/(2 + t)] leading to the optimal extracted 
work in the z-th cycle 

Wr = {htft/[2{2 + t)]. (15) 

This work still depends on the result of all measurements {-^j}i<j<i- Conditionally 
averaging this work over the last measurement will lead to a useful recursion relation as 
follows. With 

{{htf)x.^ I dX, {htfp{X,) (16) 



W. ^ hm {Wn = 7^^^ lim ({y~? - (ytf) ■ (20) 

l^OQ 2 4 + t) i-S-OO V / 



and (ED, dlO]) and (^ we get 

{{btr)x. = (^^(^ti)' + ivif - {yff (17) 

and hence 

<^''- "'>-• = (2^'^" + 2(2T0 (<^- >' - <^-*'0 ■ ('«' 

Since the last term is independent of the outcomes of measurements, subsequent 
averaging over all previous measurements {^j}i<j<i-i (indicated by an unconstrained 
bracket (...)) leads to 

(t^D = p^('^n') + spVt) ('^- )' - '^^ ■ <^«) 

Solving this recursion in the stationary limit, z — ?■ oo, we thus obtain as the average 
work per cycle 

2 + t 
2(4 + t) i- 
The last limit is easily calculated by solving the dynamics ^ for the variance as 

y\r) = l + e-'^[y\Q)-l] (21) 

and setting y'^it) = (z/j^i)^ and y^(0) = (y^) . Using f lTT]) . and identifying y^^^ with y~ 
in the limit z — )■ oo, we get in the steady state for the variance before a measurement 
the value 

{y;r ^ hm {yrf = Ul - yl + e~'\yl - 1) (22) 

+ ^{yl + 1)' - 2e-2*(y4 + 1) + e-^\yl - 1) 
with the limiting behavior 

. -^2 ^ / ym{2tY" + (1 - 2/i)t for t ^ 

^^*^ '^ I l-exp(-2t)/(l + y^) for t^oo. ^ ^ 

Likewise, the variance after a measurement becomes 

{y*?y'i 
{y*f + y? 



\2 



{ytr ^ hm {ytf = j^^^^ (24) 
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1 + y^V (1 + y. 



for t ^ 
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Finally, the average work per cycle delivered by this information machine becomes 

w. = ^^ {iy;y - (ytr) (26) 

which is our first main result, shown in Fig. [T] The power P = W^^/t becomes maximal 
if the cycle time becomes short with P ^ 1/2 ~ l/m(V2)^^^ for ^ — ^ 0. In the long 
time limit, P ~ l/[2t(l + y^)] for t — > oo. In the special case of an infinitely precise 
measurement, we get P{t, 0) = (1 - e-2*)(2 + t)/[2t(4 + t)]. 

The efficiency of this machine follows from relating the power to the rate with which 
information is acquired through the measurements. The i-th measurement yields the 
information [5] 

dxp{x\Xi) \n[p{x\Xi)/pr{x)], (27) 

which still depends on the result of all measurements {-^j}i<j<i- By using ([8]) and Qj, 
subsequent averaging over the last measurement Xi yields 

X, ^ y dX, p{X,)x, = Hyr/yt). (28) 

This simple result involves, a posteriori not surprisingly, just the variances before and 
after the measurement which are independent of the specific results {^j}i<j<j-i- ^i thus 
represents the information averaged over all measurement outcomes. In the stationary 
limit, one gets 

r (t/2)i/Vy„ - t/2 for t^O 

X. . Hy:/yt) H 1 m f i + ^ )- ^^bL^ for t -. oo. ^''^ 

I 2 V yl) 2(1+2/^)2 
Consequently, the efficiency becomes 

2(4 + t) \n{y, /yt) 
with the limiting behaviour 

f y^{t/2f/''-ylt/2 for t^O 

1-2/t (31) 

for t -> oo. 



r] 



ira) 



(l + y^)ln(l + l/2/, 

As shown in Fig. [Tl the efficiency increases monotonically with the cycle time t. It 
becomes zero for t — ;■ 0, which implies that this machine has vanishing efficiency at 
maximum power. The somewhat counterintuitive monotonic increase of 77 with the 
measurement error y^ arises from the fact that it is impossible to retrieve all information 
just by moving the center of the trap. Therefore, better measurements lead to a higher 
power but not to a higher efficiency. Indeed, while the work W^ is bounded by 1/2 
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6 



the information X^ diverges in the hmit of infinitely precise measurements y^ ~^ 0, 
leading to a vanishing efficiency. In the limit 1/^^00, both W^, and X^, tend to and 
7] — > (2 + 1)/(4 + t). For t — !■ 00, this machine can reach the upper bound 1 imposed on 
7] by thermodynamics. However, this high efficiency is somewhat useless, since in this 
case the machine delivers vanishing power. 



^ 




?s^ 




Figure 1. Performance of the machine with constant stiffness k = 1 and optimally 
controlled center A(r). Extracted work W, and efficiency 77 both as function of cycle 
time t and measurement error y^- 



For a more powerful machine, we turn to a second variant where we allow additional 
control over the stiffness of the trap A;(r). In this case, the contribution flT^ no longer 
vanishes. It becomes maximal for a standard deviation ?/(r) increasing linearly from 
y{0) = yf to 

yt + UY + i'i + tW 



y{t) = vl 



(32) 



'^+1 2 + t 

In the stationary limit, z — )■ 00, using f l32|) instead of fl2T]) and the same reasoning to 
derive the limiting behavior as above, we obtain for the variance prior to a measurement 
in the steady state {y^Y the cubic equation 



[{2 + t){y:f-ty[{y 
The limiting behavior of its solution is 

y^ei^ + {?,/A-yl)t/2 

-7- 



4(y. 



)'yl 



0. 



{yl 



for t ^ 

for t — 7- 00. 



(33) 



(34) 



{^ + ylY'\ 

For ( !24l) one obtains {ytY with the short time and quasistatic behavior 

yrJ^/^-{b/A + yl)t/2 fort 



{yt 



+ \2 



%m(l+^m-l/m(l+ y; 

(1 + ylYt 



2 Nl/2^ 



for t — !■ 00. 



(35) 
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For this second variant, we can still determine the contribution to the extracted work 
from flT3|) as in the first case, provided we use the solution of fl33|) in the expression fl20l) 
for the stationary limit. Collecting everything, we obtain for the extracted work the 
expression 

^^ ^ _hh?__{gf_ _ {j£_yt)_ ^ in(^-/^+) (36) 

shown in Fig. |2l In this case, the power diverges in the short limit as 

whereas in the long time limit one obtains 

P«iln(lH-±). (38) 

The short time divergence of the power is compensated by a corresponding divergence of 
the rate of information acquired through the measurements. Indeed, similarly as above, 
one gets the information per measurement 

r t^/V(22/™) - (4 + l/yl)t/lQ for t ^ 

x.^Hy:lyt)^\ li,A^M_i-W(i + .^)^ for t ^ J.^^^ 
I 2 V yl) t{^ + yl) 

The efficiency of this machine becomes 

{{y:f-{yt?)l{^ + t) + {y:-ytfit 

\n{y-/yt) 
t'/yi8ym) for t^O 




(40) 



4(1-W(1 + .^)V-) ^^^ ^^ (41) 

tHl + l/yl) 

shown in Fig. |2J 

For this variant, the efficiency increases with the cycle time starting at r^ = 1/2 for 
t — > and saturating the upper bound 77 = 1 for t — )■ oo and any ym- In this quasistatic 
case, in contrast to the first variant, the two control parameters allow to extract the 
full information. In another difference, the efficiency monotonically decreases with 
increasing ym- Here, more precise measurements lead to a larger efficiency allowing 
even 77 = 1 at finite t for infinite precision ym — ;• 0. In the full (t, |/m)-plane, the 
efficiency is bounded by 1/2 from below. The value 1/2 found here in the short time 
limit that corresponds to maximum power may hint to a relation of our result with 
that for the efficiency of isothermal machines at maximum power where the value 1/2 
is universal in the linear response regime [IT]. While it is not obvious how to map 
repeated measurements for short cycle times to a linear response formalism, finding the 
same value in both cases may be more than incidental. 

In conclusion, we have studied the efficiency for a cyclically operating Brownian 
information machine consisting of an overdamped particle in a time-dependent harmonic 
trap. For two variants of such a machine, we have obtained analytically how the 
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Figure 2. Performance of the machine with optimally controlled stiffness fc(r) and 
center A(t). Extracted work Wtt and efficiency r\ both as function of cycle time t and 
measurement error y„i. 



efficiency depends on both the precision of a positional measurement and the cycle 
time. Beyond these specific results our work raises a few questions concerning such 
machines in general. First, while the quite natural definition of efficiency defined as 
mean extracted work divided by the mean acquired information shares features such as 
boundedness between and 1 with the more conventional thermodynamic definition of 
efficiency for ordinary isothermal machines, finding 77 = 1 even for finite cycle time in 
the limit of infinitely precise measurements, as we do for the second variant, suggests 
that these information machines differ in essential aspects from thermodynamics ones. 
For reaching ?7 = 1, the latter require a quasistatic operation, i.e., an infinite cycle time. 
Second, is it possible to formulate a linear response theory, i.e., to calculate Onsager 
coefficients for such machines? Third, can we derive general bounds on the efficiency 
at maximum power following reasoning for non-feedback driven machines? Finally, an 
experimental test of such a machine would be interesting and should be possible with 
available technology. 
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